Constraints on Neutrino-Nucleon Interactions at energies of 1 EeV with the IceCube 

Neutrino Observatory 
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A search for extremely high energy cosmic neutrinos has been carried out with the IceCube 
Neutrino Observatory. The main signals in the search are neutrino-induced energetic charged leptons 
and their rate depends on the neutrino-nucleon cross section. The upper-limit on the neutrino flux 
has implications for possible new physics beyond the standard model such as the extra space-time 
dimension scenarios which lead to a cross section much higher than the standard particle physics 
prediction. In this study we constrain the neutrino-nucleon cross section at energies beyond 10 9 GeV 
with the IceCube observation. The constraints are obtained as a function of the extraterrestrial 
neutrino flux in the relevant energy range, which accounts for the astrophysical uncertainty of 
neutrino production models. 

PACS numbers: 98.70.Sa, 95.85.Ry, 13.15.+g 



I. INTRODUCTION 

High energy cosmic neutrino observations provide a 
rare opportunity to explore the neutrino-nucleon (i>N) 
interaction behavior beyond energies accessible by the 
present accelerators. These neutrinos interact during 
their propagation in the earth and produce energetic 
muons and taus. These secondary leptons reach under- 
ground neutrino detectors and leave detectable signals. 
The detection rate is, therefore, sensitive to neutrino- 
nucleon interaction probability. The center-of-mass en- 
ergy of the collision, y 7 ?, is well above ~ 10 TeV for 
cosmic neutrino energies on the order of 1 EeV (= 10 9 
GeV). This is a representative energy range for the bulk 
of the GZK cosmogenic neutrinos, generated by the inter- 
actions between the highest energy cosmic ray nucleons 
and the cosmic microwave background photons [l|. 

The i/N collision cross section can vary greatly if non- 
standard particle physics beyond the Standard Model 
(SM) is considered in the high energy regime of y/S 
TeV. The extra-dimension scenarios, for example, have 
predicted such effects 0, In these scenarios, the vir- 
tual exchange of Kaluza-Klein graviton Q, or micro- 
scopic black hole production Q leads to a substantial in- 
crease of the neutrino-nucleon cross section by more than 
two orders of magnitude above the SM prediction. The 
effect would be sizable enough to affect the expected an- 
nual event rate (0(0.1 — 1)) of the GZK neutrinos in the ~ 
km 3 instrumentation volume of an underground neutrino 
telescope such as the IceCube observatory. Thereby the 
search for extremely-high energy (EHE) cosmic neutrinos 
leads to constraints on non-standard particle physics Q . 

The IceCube neutrino observatory has already begun 
EHE neutrino hunting with the partially deployed un- 
derground optical sensor array @. The 2007 partial Ice- 
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Cube detector realized a ~ 0.7 km 2 effective area for 
muons with 10 9 GeV and recently placed a limit on the 
flux of EHE neutrinos approximately an order of mag- 
nitude higher than the expected GZK cosmogenic neu- 
trino intensities with 242 days of observation [7J. Since 
new particle physics may vary the cross section by more 
than an order of magnitude as we noted above, this re- 
sult should already imply a meaningful bound on the i/N 
cross section. In this paper, we study the constraint on 
the cross section (ct^n) by the null detection of EHE 
neutrinos with the 2007 IceCube observation. A model- 
independent bound is derived by estimating the lepton 
intensity at the IceCube depth with the SM cross section 
scaled by a constant. The constraint is displayed in the 
form of the excluded region on the plane of the cosmic 
neutrino flux and <t„n . It is equivalent to an upper-bound 
on er„N for a given flux of astrophysical EHE neutrinos. 
We also study the mo del- dependent constraint on the 
microscopic black hole creation by neutrino-nucleon col- 
lision predicted in the extra-dimension scenario @. We 
calculate the fluxes of leptons propagating in the earth 
including the black hole cross section and the final states 
to estimate expected event rate in an equivalent IceCube 
2007 measurement as a function of extraterrestrial neu- 
trino intensity. The null detection of signal candidates 
leads to a constraint on this particular scenario. 

There are several works on model-independent upper 
bounds of ovn using the observational limit of EHE neu- 
trino flux in the literature. Refs. [!, 0, [|| derived the 
bound using the results of horizontal air shower search 
by AGASA [| and Fly's Eye [l0|. Refs. d B3 set the 
limit based upon the flux bound by the RICE experi- 
ment [lH • Our approach in the present study is different 
mainly in two respects. The previous works assumed the 
GZK cosmogenic neutrino bulk as the guaranteed beam 
and deduced the cross section limit using the GZK neu- 
trino intensity. Here extraterrestrial neutrino intensity is 
considered as a free parameter. This method is an ap- 
plication of the technique to derive the flux limit based 
upon the quasi-differential event rate @, [IH, [l3[ , which 



2 



is independent of specific neutrino flux models. As EHE 
cosmic ray composition and their origin are still quite 
uncertain, this approach provides more appropriate con- 
servative limits on ct^n ■ It also allows estimation of the 
minimum intensity of neutrino flux required to constrain 
the cross section. Another difference is that the previ- 
ous works introduced the simplification that event rate 
solely depends on rates of electromagnetic or hadronic 
cascades directly initiated by neutrinos inside the ef- 
fective volume of the detector. This is in fact a good 
approximation for the RICE experiment which is sen- 
sitive to radio emission from shower events. However, 
underground neutrino telescopes such as IceCube have 
larger effective area for through-going muons and taus in 
EHE neutrino search [T3 - [l6{ . This study of the model- 
independent limit includes calculation of not just inten- 
sities of neutrinos but also the secondary muon and tau 
fluxes reaching the detection volume for a given ovn and 
includes their contributions in the overall event rate. 

The paper is outlined as follows: First we discuss the 
model-independent constraint in Sec. |TTJ The method to 
calculate the neutrino and the secondary lepton propa- 
gation from the earth's surface to the IceCube detector 
depth is described. The fluxes for different strengths of 
oVn are calculated and the resultant constraint is shown 
for both cr„N and the cosmic neutrino flux at neutrino en- 
ergies of 1 and 10 EeV, respectively. Sec. lIIII describes the 
constraint on the microscopic black hole production by 
neutrino-nucleon interaction as an example of the model 
dependent bound on ovn- Fluxes of muons and taus 
from evaporation of black holes produced in the neutrino- 
nucleon collision in the earth are calculated. Their contri- 
butions, as well as those from contained hadronic showers 
induced directly by the evaporation, would give an ob- 
servable event rate in the IceCube 2007 measurement, 
and thereby put constraints on the black hole scenario. 
We summarize our conclusions in Sec. IIVI 



II. MODEL INDEPENDENT CONSTRAINT ON 
THE NEUTRINO-NUCLEON CROSS SECTION 

The flux limit obtained by the present IceCube ob- 
servations allows us to place an upper bound on the 
neutrino-nucleon cross section in a model independent 
manner; new physics cannot increase ovn too much, oth- 
erwise EHE neutrinos would have produced observable 
events. As an underground neutrino telescope is sensi- 
tive to not just shower events induced from neutrinos, 
but to through-going muons and taus generated by the 
neutrino-nucleon scattering, one must understand how 
much fluxes of these leptons reaching an underground 
detection volume is increased with ovn- In this section, 
we first discuss our method to calculate intensities of neu- 
trinos, muons, and taus at the underground depth of the 
IceCube observatory for a wide range of ovn strength, 
followed by a description of how they would contribute 
to the event rate. Finally the constraint on both a„N 



and cosmic neutrino flux is described together with the 
relevant discussions. 



A. The Method 

Given a neutrino flux at the surface of the Earth, 
the neutrino and charged lepton fluxes at the IceCube 
depth are calculated by the coupled transportation equa- 
tions fig : 
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where Ji = dNi/dEi and J v = dN v /dE u are differential 
fluxes of charged leptons (muons and taus) and neutri- 
nos, respectively. X is the column density, Na is the 
Avogadro's number, p is the local density of the medium 
(rock/ice) in the propagation path, a is the relevant inter- 
action cross section, dnf jdE is the energy distribution of 
the decay products which is derived from the decay rate 
per unit energy, c is the speed of light, m/ and t, are the 
mass and the decay life time of the lepton I, respectively. 
CC(NC) denotes the charged (neutral) current interac- 
tion. In this study we scale cj„n to that of the SM pre- 
diction with the factor N sca i c , i.e., ct„n = N sca i e cr^. It 
is an extremely intensive computational task to resolve 
the coupled questions above for every possible value of 
ovn. To avoid this difficulty, we introduce two assump- 
tions to decouple calculation of J v from the charged lep- 
ton transportation equation. The first is that distortion 
of the neutrino spectrum by the neutral current reaction 
is small and the other is that regeneration of neutrinos 
due to muon and tau decay and their weak interactions is 
negligible. These are very good approximations in the en- 
ergy region above 10 8 GeV where even tau is unlikely to 
decay before reaching the IceCube instrumentation vol- 
ume. Then the neutrino flux is simply given by the beam 
dumping factor as 

J^E^Xjc) = J w (E„,0)e- JV - 1 ^N' oo Xi C) ( 3) 

where Xic is column density of the propagation path 
from the earth surface to the IceCube depth. The 
charged lepton fluxes, J; =/1)T (Ei, Xic), are obtained as 

J M , T (E M , r ,X IC ) = N A f dX j dE^ T ^- (F4 jT E M , T ) 
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FIG. 1: Integral fluxes of the muon and taus above 10 PeV 
(= 10 7 GeV) at IceCube depth (~ 1450m) for GZK cosmo- 
genic neutrinos [17[ ■ The solid lines represent muons while the 
dashed lines represent taus. Numbers on each of the curves 
are the multiplication factors (N aca i e ) that enhance the stan- 
dard fN cross section flol in the relevant calculations. 
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Here dN^^ T /dE tl ^ T (F lfi T — !• E^.i-) represents distributions 
of muons and taus with energy of E M!r at Xic created by 
2/N collisions at depth X with an energy E p T . This is 
calculated in the transportation equation, Eq. [21 with a 
replacement of J U (E V ) by Eq. [3J 

Calculation of the neutrino and the charged lepton 
fluxes with this method is feasible for a wide range of 
Nscaic without any intensive computation. A compari- 
son of the calculated fluxes with those obtained without 
the introduced simplification for a limited range of N sca j e 
indicates that the relative difference we found in the re- 
sultant Jv,il,t(Xjc) is within 40%. Since this analysis 
involves an order of magnitude of increase in er„N , the in- 
troduced approximations provide sufficient accuracy for 
the present study. 

FIG.[TJshows the calculated intensities of the secondary 
muons and taus for various N sca i e factors. Here the pri- 
mary neutrino spectrum is assumed to follow the GZK 
cosmogenic spectrum and flux calculated in Ref. [l7[ as- 
suming an all-proton cosmic-ray composition with a mod- 
erately strong source evolution, (1 + z) m with m = 4 
extending to z = 4. One can see that the intensity is 
nearly proportional to N SC aic as expected since the inter- 
action probability to generate muons and taus linearly 
depends on ovn- It should be pointed out, however, that 
the dependence starts to deviate from the complete lin- 
earity when the propagation distance is comparable to 
the mean free path of neutrinos, as one can find in the 
case of N sca io = 10 in the figure. This is because the neu- 
trino beam dumping factor in Eq. [3] becomes significant 
under this circumstances. 

The flux yield of leptons at the IceCube depth, Yl 
(I = v s, /i,T), originating from neutrinos with the same 
energy at the earth's surface, Ef,, is given by Eq. @]for 



muons and taus, by Eq.|3]for neutrinos, with an insertion 
of J V {E U , 0) = 5{E V — El). The resultant event rate per 
neutrino energy decade is then obtained by @, [H, G3 > 



dJ Vc 



3 dlogEv 



-{Et 



dQ 



£ JdEiAi{Et) Y^(El,Ei,Xic(ii),N s 



a) (5) 



where A\ is the effective area of the IceCube to detect the 
lepton I. In this equation above, the I = fj,,T terms rep- 
resent the through-going track events while contribution 
of events directly induced by neutrinos inside the detec- 
tion volume is represented by the terms I = u ei z^, v r . 
The effective area for v s, A Ul is proportional to ovn i.e., 
Nscaio so the rate of contained shower events is linearly 
dependent on the neutrino-nucleon scattering probabil- 
ity. Note that the differential limit of the neutrino flux 
is given by Eq. [5] for N sca i c — 1 with N u = jlgo which 
corresponds to the 90 % confidence level average upper 
limit. This calculation is valid when the cosmic neutrino 
flux J v and the cross section ovn do not rapidly change 
over a decade of neutrino energy around Ef,. Limiting 
ovn in the present analysis corresponds to an extraction 
of the relation between N sca i c and the (unknown) cos- 
mic neutrino flux J v ^+ v ^+ Vt yielding N v = /2go- The 
obtained constraints on ct„n is represented as a function 
of J Vc +v +v T for a given energy of E s v . It consequently ac- 
counts for astrophysical uncertainties on the cosmic neu- 
trino flux. 

In scenarios with extra dimensions and strong gravity, 
Kaluza-Klein gravitons can change only the neutral cur- 
rent (NC) cross section because gravitons are electrically 
neutral. Any scenarios belonging to this category can be 
investigated by scaling only in the present analysis. 
The event rate calculation by Eq.[5]is then performed for 
Vj (N sca i c = 1) with effective area for z/'s, A v , enhanced 
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the rate of detectable events via the NC reaction by Ice- 
Cube is proportional to . We also show the con- 
straint in this case. 



B. Results 

In this analysis we use the IceCube observation results 
with 242 days data in 2007 to limit ovn using Eq. [5] No 
detection of signal candidates in the measurement has led 
to an upper limit of the neutrino flux of 1.4 x 10~ 6 GeV 
cm~ 2 sec -1 sr" 1 7] in the energy range from 3 x 10 7 to 
3 x 10 9 GeV. The effective area A t is ~0.7 km 2 for fi, ~ 
0.4 km 2 for t, and 3x 10~ 4 km 2 for v s p}. Constraints on 
ovn are then derived with Eq.[5] The results for E s v = 10 9 
and 10 10 GeV are shown in FIG. [2] Enhancing the 
charged current cross section by more than a factor of 30 
for E v = 1 EeV (10 9 GeV) is disfavored if the astrophys- 
ical neutrino intensities are around ~ 10~ 7 GeV cm~ 2 
sec -1 sr -1 , near the upper bound of the GZK cosmogenic 
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FIG. 2: Constraints on the all-flavor sum of cosmic neutrino 
flux and the charged current vN cross section based on the null 
detection of neutrino signals by the IceCube 2007 observation. 
The right upper region is excluded by the present analysis. 
The cross points provide reference points where the standard 
cross section and the expected GZK cosmogenic neutrino 
fluxes [3 is located. 
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FIG. 3: Constraints on the all-flavor sum of cosmic neutrino 
flux and the neutral current i^N cross section for the scenario 
that only the neutral current reaction is enhanced by a new 
physics beyond the standard model. The right upper region 
is excluded by the present analysis. The crosses provide ref- 
erence points for the standard cross section [Ts| and the ex- 
pected GZK cosmogenic neutrino fluxes fisl j - 



neutrino bulk. Note that neutrino-nucleon collision with 
E u = 1 EeV corresponds to yfs ~ 40 TeV and the present 
limit on <7„n would place a rather strong constraint on 
scenarios with extra dimensions and strong gravity, al- 
though more accurate estimation requires studies with a 
mo del- dependent approach which implements the cross 
section and the final-state particles from the collision pre- 
dicted by a given particle physics model. Taking into 
account uncertainty on the astrophysical neutrino fluxes, 
any model that increases the neutrino-nucleon cross sec- 
tion to produce charged leptons by more than two orders 
of magnitude at y/s ~ 40 TeV is disfavored by the Ice- 
Cube observation. However, we should point out that 
the IceCube 2007 data could not constrain the charged 
current cross section if the intensity of cosmic neutrinos 
in the relevant energy region is fewer than ~ 10~ 8 GeV 
cm sec -1 sr , within the lower range of prediction for 
the cosmogenic neutrino fluxes [2(j. Absorption effects 
in the earth becomes sizable in this case, resulting in less 
sensitivity to the cross section. This limitation will be 
improved for larger detection area of the full IceCube 
detector. 

FIG. [3] shows the constraints when only the NC cross 
section is varied. Enhancement of by a factor be- 
yond 100 at ~ 40 TeV is disfavored, but this strongly 
depends on the cosmic neutrino flux one assumes. Be- 
cause the NC interaction does not absorb neutrinos dur- 
ing their propagation though the earth, the cross section 
could be bounded even in the case when the neutrino flux 
is small, but the limit becomes rather weak; the allowed 
maximum enhancement factor is on the order of ~ 10 3 . 



III. CONSTRAINT ON THE MICROSCOPIC 
BLACK HOLE PRODUCTION 

A constraint on a specific physics model that enhances 
the neutrino-nucleon cross section is obtained by the 
same procedure for the model independent bound, ex- 
cept the transport equations, Eqs. Q] and [U would have 
total and differential neutrino cross section provided by 
both SM and the new model. Here we study the model of 
black hole creation as a possible consequence of low-scale 
gravity that may occur if space-time has more than four 
dimensions. We use the predicted cross section of lack 
hole production via the neutrino-nucleon scattering de- 
scribed by Ref. @, parametrized by the Planck scale Mrj, 
the ratio of the minimal black hole mass to the Planck 
scale x m i ni and the space-time dimension D = 4 + n. 
In this paper n = 6, and Mjj = ITeV are assumed as 
representative numbers. The resultant cross section may 
exceed SM interaction rates by two orders of magnitude 
or even greater. Therefore, the model independent bound 
shown in the previous section indicated that the 2007 Ice- 
Cube observation should be already sensitive to some of 
the parameter space in the black hole creation model. 

The final states in the neutrino-nucleon scattering in 
this model are quite different from the SM case. Black 
holes evaporate and generate multiple particles of all 
kinds, like leptons, quarks, gluons, and bosons. These 
products are distributed according to the number of de- 
grees of freedom. Consequently, the average number of 
muons and taus, N^ +T , are 1/30 of all particle average 
multiplicity N, which is also determined by the specific 
model. As N ~ 10 at neutrino energy of E v = 1 EeV, 
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FIG. 4: Left: Integral fluxes of the muon and taus above 10 PeV (= 10 7 GeV) at IceCube depth (~ 1450m) of the GZK 
cosmogenic neutrinos [13] in case of the microscopic black hole creation scenario [fj for (Mj,a; m i n ) = (lTeV,l) (solid), (lTeV,3) 
(dot-dash). The dash line corresponds to the intensities obtained by the SM uN cross section [lq]. Right: Energy spectra of 
the GZK v induced muons and taus at IceCube depth with downgoing (i.e. cos(zenith) > 0) geometry expected by the black 
hole model for (Mo,a; TO j rl )=(lTeV,l). The spectra produced by the SM cross section are also shown as the dashed lines for 
comparison. The curve labeled "Atm /x" represents the atmospheric muon intensity estimated by the IceCube observation with 
its uncertainty expressed by the two dash lines [2l| . 



multiple muon or tau production would very rarely occur. 
Then the effective differential cross section da u a/dE^ r 
in the transport equation Eqs. [T] and [5] in the black hole 
model is represented by 

dEZ T = 2 U ^ {Ev) ^lEV (6) 

with < E^ T < 2E V /N. We take N from Ref. [5] in the 
present calculation. In this specific scenario, a muon or 
a tau carries a small fraction (1/N ~ 0.1) of incoming 
neutrino energy E v in average, in contrast to the SM 
collision that takes away 1 — y ~ 0.8 of neutrino energy 
by a generated charged lepton. 

Solving the transport equations gives the intensities of 
secondary muons and taus, which are shown in FIG. 2] 
One can find in the zenith angle distribution (the left 
panel) that the intensities are increased by more than 
two orders of magnitude above the SM case. The large in- 
crease of er„N enhances down-going event rates while the 
upgoing muon and tau rates are more suppressed. The 
zenith angle distribution is consistent with the original 
work in Ref. It should also be noted that the energy 
spectra is substantially modified from those in the SM 
case (the right panel). The peak energy is around 1 EeV, 
an order of magnitude lower than the SM spectrum, re- 
flecting the fact that a smaller fraction of neutrino energy 
is channeled into muons and taus via the black hole evap- 
oration. The peak happens to match the most sensitive 
energy region in the IceCube EHE neutrino search Q- 

Because ovn is solely predicted by the specific model, 
the model-dependent constraints on z/N interactions is 
represented in the plane of extraterrestrial neutrino flux 
and the number of events the IceCube 2007 run would 
have detected. FIG. [5] shows the number of events as 




Number Of Events 

FIG. 5: Number of events as a function of extraterrestrial neu- 
trino flux at 1 EeV for (M D ,x min )=(lTeV,l) (solid), (lTeV,3) 
(dash). The 90 % C.L. line determined by the Poisson statis- 
tics is also shown as the vertical line for reference. 



a function of the neutrino intensity at energy of 1 EeV, 
if the microscopic black hole evaporation occurs as in 
Ref. Q. The Poisson statistics then determines the 
upper-limit of neutrino flux that can be still consistent 
with the null observation by IceCube. It is indicated that 
the neutrino intensity of 10 -7 GeV cm~ 2 sec -1 sr _1 is 
disfavored in this scenario. More parameter space of Md 
and x m in will be further constrained by near future ob- 
servation with IceCube whose detection volume is rapidly 
growing with increase of number of the detectors in op- 
eration. 
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IV. CONCLUSIONS 

The IceCube 2007 observation indicated that any sce- 
nario to enhance either the NC or both the NC and CC 
equivalent cross section by more than 100 at y/s ~ 40 
TeV is unlikely if sum of the all three flavors of astro- 
physical neutrino fluxes are greater than ~ 3 x 10~ 8 
GeV cm -2 sec -1 sr _1 in EeV region. Many models of 
the GZK cosmogenic neutrinos exist to predict this flux 
range, thus the present constraints limit new particle 
physics beyond SM, unless the extraterrestrial neutrino 
intensity is smaller than expectation. The example of 
the mo del- dependent bound on ctvn has been also shown 
for the microscopic black hole evaporation scenario. A 
high cosmic neutrino intensity constrains the parame- 
ter space of the black hole creation. Future observation 
by the rapidly growing IceCube detectors will strongly 



limit particle physics models which predict an increase 
of neutrino-nucleon interaction probability. 
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